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The gravitational action of the smooth energy-matter components filling in the universe can affect the orbit 
of a planetary system. Changes are related to the acceleration of the cosmological scale size R. In a universe 
with significant dark matter, a gravitational system expands or contracts according to the amount and equation 
of state of the dark energy. At present time, the Solar system, according to the LCDM scenario emerging from 
observational cosmology, should be expanding if we consider only the effect of the cosmological background. 
Its fate is determined by the equation of state of the dark energy alone. The mean motion and periastron 
precession of a planet are directly sensitive to R/R, whereas variations with time in the semi-major axis and 
eccentricity are related to its time variation. Actual bounds on the cosmological deceleration parameters go from 
accurate astrometric data of perihelion precession and changes in the third Kepler's law in the Solar system fall 
short of ten orders of magnitude with respect to estimates from observational cosmology. Future radio-ranging 
measurements of outer planets could improve actual bounds by five order of magnitude. 
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I. INTRODUCTION 



The local effect on particle dynamics from the global ex- 
pansion of the universe is a classical topic in general rela- 
tivity. The picture is pretty clear for a test particle moving 
in the homogeneous and isotropic Robertson- Walker (RW) 
metric. Analytical methods are mainly based on solving the 
geodesic equations JH and references therein] and have been 
used to examine whether particles removed from the Hubble 
flow asymptotically rejoins it. Different definitions of what 
this means can lead to apparent disagreement [2]. A general 
result for a test particle in a RW metric is that its peculiar ve- 
locity with respect to the Hubble flow declines in magnitude as 
the scale factor of the universe, R, expands, alike to any pho- 
ton energy [3]. This characteristic does not assure that all of 
the features of a generic trajectory in the expanding space will 
approach those of a trajectory in the Hubble flow. As a matter 
of fact, many characteristics depend on the acceleration, or de- 
celeration, of the universe rather than just its expansion. The 
strong requirement that the proper distance between the test 
particle and a reference asymptotic particle approaches zero 
holds only in an accelerated universe. In an Einstein-de Sitter 
universe, a particle initially at rest with respect to the origin 
falls towards the origin, passes through it and asymptotically 
moves on the opposite side of the sky with a null peculiar ve- 
locity JUS]- Anyhow, the spatial 'peculiar' distance between 
the test particle and the background particle with the same 
asymptotic speed is unbounded and increases with time fl. 
Even for an under-dense universe, the particle always keeps 
a distance from its corresponding background place 0]. A 
cosmological constant can change this picture. Due to the re- 
pulsive effect of the vacuum energy, the peculiar distance de- 
creases asymptotically. Nevertheless, the effect does not arise 
from the expansion of the universe but stems from the gravi- 
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tational repulsion of the new ingredient in the cosmic energy 
budget Jill. 

Matter inhomogeneities complicate the picture. In the 
Einstein-Strauss model, the strict vacuum spherically sym- 
metric Schwarzschild solution is truncated at a finite radius 
and matched to a cosmological model. The Swiss cheese 
universe generalizes this solution, prescribing that the mass 
within non overlapping spherical cavities in an otherwise ho- 
mogeneous and isotropic expanding universe is compressed 
to a point mass. During the expansion the only change of any 
cell is an increase of the truncation radius but, within the cav- 
ity, the gravitational field remains exact Schwarzschild and the 
behavior of the rest of the universe is irrelevant. This could 
be a compelling argument on the (absence of) cosmological 
influence on gravitationally bound systems, its obvious short- 
coming being that it does not consider the gravitational ef- 
fect of the (expanding) energy-matter fluid wherein the test 
particles move. In analogy with the RW metric, it should be 
not the expansion of the universe that can cause gravitational 
bound systems to separate or evolve. The effect, if any, should 
stem from the gravitational action of the homogeneous back- 
ground in which the system is embedded. How the expansion 
of the universe affects the dynamics of local systems has been 
analyzed mainly considering the two-body problem in a cos- 
mological background 

H IE H • Some apparent differ- 
ing results could spring from the different coordinate systems 
adopted in different studies. The equation of motion are not 
coordinate invariant and the form of any correction due to the 
cosmic background is dependent upon the frame employed 
01 • I n f act ' many apparent discrepancies disappear if the dif- 
ferent approaches are performed coherently [10]. 

In a series of previous papers, we discussed the effect on 
a local bound system of the cosmological constant ifTll [l2ll 
and of diffuse dark matter 11311 . In this paper we investigate 
the evolution of a planet orbit around a massive star embed- 
ded in an homogeneous and isotropic universe. In Section fill 
we use some arguments based on Newtonian dynamics to re- 
derive the relevant equations of motion for a test particle in a 
RW universe perturbed by a weak local inhomogeneity. Sec- 
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tion[TTT]is devoted to study the effect of dark matter and dark 
energy on the orbital radius evolution of a bound system. In 
Section HVl the Lagrange's planetary equations are discussed 
and the time variation of the orbital elements is derived. Sec- 
tion [V] is devoted to some final comments. 



II. EQUATION OF MOTION 

The essential physics of the motion of a test particle in the 
weak gravitational field of a local inhomogeneity perturbing 
an otherwise homogeneous and isotropic expanding universe 
can be extracted following the Newtonian approach as usually 
performed for the dynamics of linear perturbations flill . The 
background is suitably described by the RW metric. The well 
known equations for the evolution of the scale factor R are 



Rr = 
R = 



8irG 



, P R 2 -kc 2 , 
—-R{p + 3—), 



(1) 
(2) 



where dot stands for time derivative, G is the Newton's grav- 
itational constant, c the speed of light in the vacuum, p and p 
the density and pressure of the total cosmic fluid (inclusive of 
any dark energy component), respectively, and k a curvature 
constant. Equation (Q~|) is known as the Friedmann equation 
and can be also seen as an equation for the Hubble parameter 
H = R/R. The scale factor R can be expressed in terms of 
the redshift z as R = 1/(1 + z), where we have normalized 
R = 1 at the present time. 

Let us introduce Eulerian coordinates x. The corresponding 
comoving spatial coordinates are r = x(t)/R(t). The time 
derivative of x is the sum of the Hubble expansion plus the 
peculiar velocity, x = Hx + <5v, where the peculiar velocity 
<5v can be expressed as the time derivative of the comoving 
coordinate, 5v = R r. Differentiating x twice gives 



R 

x = <5v + H5v H x. 

R 



The Newtonian limit of the equation of motion for a test parti- 
cle in the RW metric is obtained from Eq. 01 for zero peculiar 
velocity and zero peculiar acceleration, x = (R/R)x. 

The perturbed equation of motion follows from the Eulerian 
equation of motion 



go + g, 



(4) 



with go being the unperturbed gravitational acceleration that 
acts on a particle in a homogeneous universe and g the pecu- 
liar gravitational acceleration. We remind that go can not be 
derived from Newtonian gravity alone and must be assumed 
to be given IU411 . 

Identifying the right hand sides of Eqs. in the case 
of zero peculiar terms, we get the unperturbed equation, 
(R/R)x — go. Then, subtracting such an unperturbed equa- 
tion from Eq. (0 side by side, we get 



R 

7r 



g- 



(5) 



The cosmological influence on the local dynamics of a bound 
system appears through the acceleration of the scale factor and 
not through its first derivative. Equation 01 has been already 
derived several times in literature with g specified to the case 
of a point mass, g = — GMx/x 2 151 1 1 Oil . Cooperstock et al. 
J5t] derived the geodesic equations for a RW metric in the lo- 
cal inertial frame using Fermi normal coordinates and then 
added linearly the peculiar Newtonian acceleration due to a 
massive body. Carrera and Giulini [ 10] reviewed different ap- 
proaches leading to Eq. (|5j. In particular they considered both 
the motion of a test particle in the McVittie spacetime, al- 
ready considered in 11511 . and the full relativistic treatment for 
electromagnetically-bounded systems, drawing on some argu- 
ments first discussed in ]7|]. 

With the aid of Eq. 0), equation 0) can be rewritten in 
terms of the peculiar velocity as 



5v + H5v = g. 



The solution of equation (0 is 



5v(t) 



R(t) 



A 



R(t')g(t')dt' 



(6) 



(7) 



with A a constant of integration. In the unperturbed case, pe- 
culiar velocities of non-relativistic objects declines as 1/ R(t), 
just as any photon momenta. Of course, in presence of a local 
inhomogeneity in the matter distribution, Sv does not tend to 
zero anymore. Dot in the previous equations stands for the 
full convective time derivative, which when dealing with per- 
turbed quantities can be replaced by d/dt = d / d t + vq ■ V. 
For a point-mass perturbation in a uniform background, 
d(8po/po)/dt = 0, and the perturbed linearized equation for 
conservation of energy takes the form V<5v = 0. Then, the 
convective derivative can be identified with the partial time 
derivative, d5v/dt = dSv/dt. 

In the following sections, Equation (0 will be specified in 
the case of a central massive body, g = —GMic/x 2 , and will 



(3) be the starting point for our considerations. 



III. EVOLUTION OF THE ORBITAL RADIUS 

The motion of a test particle in the expanding universe is a 
typical multiple time scale problem in which the natural pe- 
riod of the unperturbed orbit is much smaller than the time 
scale over which the perturbation changes, i.e. the age of the 
universe. Equation (0 can be conveniently rewritten in cylin- 
drical coordinates {x, ip} JUdll]. For a test particle around a 
massive central object M, it is 



xip 



GM R 
— t, -x = 



R 

x 2 tp 



io- 



CS) 
(9) 



with Lq the conserved angular momentum per unit mass. 
Since the perturbation is purely radial, the motion is confined 
to the unperturbed Keplerian orbital plane. An approximate 
solution to the orbital motion can be found if we consider the 
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FIG. 1: Orbital radius evolution for an Earth-like planet (orbital pe- 
riod ~ 1 year). The variation is put to zero at the transition redshift 
z t i from a decelerated to an accelerated expansion. We consider flat 
models with h = 0.7. The full line is for an Einstein-de Sitter model 
Omo = 1 (the acceleration gets asymptotically to zero at zti — ► —1), 
the dashed line for a LCDM model with Qmo ~ 0.3 and wxo = — 1 
(Ztr — 0.67), the long-dashed line for Qmo = 0.3 and phantom 
energy with wxo = — 2 (z tr ~ 0.51). 
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FIG. 2: The curves represent the loci of points in the space of dark 
energy parameters, i.e. the Qxo — ffixo plane, for which the radius 
of a gravitationally bound system does not vary at the present time. 
Curves from the left to the right are for Qmo = 0.2, 0.3 and 0.4, 
respectively. The region on the left side of each curve corresponds to 
positive expansion for the given Qmo ■ 



perturbations of the orbital coordinates in the case of circular 
orbits, 



x{t) = .T + Sx(t), 
<p{t) = n Q t + 5<p(t), 



(10) 
(11) 



with xo and no — \J GMjx\ being the unperturbed orbital 
radius and mean motion, respectively. We assume that the 
time dependence of the perturbations in Eqs. dTOllTTb is due to 
the time variation of the perturbing function. Substituting in 
Eqs. dUU) yields 



c" 2 r R 

ox + n Q dx — —xq 
R 



0. 



(12) 



FIG. 3: Time (in units of 10 9 /h years) needed for the dissociation 
of an Earth-like system (P ~ 1 year) in a phantom cosmology as a 
function of the equation of state of the dark energy. We have consid- 
ered a flat model with £Imo = 0.3. 



As far as we consider secular variations, the Sx term in the 
above equation can be neglected since the unperturbed orbital 
period is much smaller than the universe age [5]. Then, 



x(t) ~ Xq I 1 + 



1 R. . 

^!« (t) , 



(13) 



This can be also seen by writing down the equation for the 
centripetal acceleration in a circular motion, 



ip 2 x 



GM R 



R 



(14) 



where <p is the angular frequency that we can assume as con- 
stant during a planet orbit. Remembering that (fix 2 = Lq — 



x 

XQ 



1 R 
~njR 



(15) 



Inserting an expansion for x such as in Eq. ( TTOb in the above 
equation, we re-obtain the result for the secular evolution, 
Eq. ( fT3l i. The unperturbed orbital radius can be fixed to its 
value at the time when the perturbation cancels out, i.e. when 
the expansion of the universe changes from a decelerated ex- 
pansion to an accelerated one. 

In what follows, we will consider a model of universe filled 
in with pressureless dark-matter, 57mo, and dark energy f2xo 
with a constant equation of state, wxo- Then the field equa- 
tions for the scale factor, Eqs. ([Tl|2|, take the form 



R 
R 



Hi 



(17) 



x {n M0 R- 3 + (1 + 3u. XO )^xoir 3(1+ ™ xo) } 



with Qko = 1 — Omo — ^xo- The Hubble constant can be 
written in convenient units as Hq = h x 100 km s _1 Mpc _1 . 
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For our reference model, we will mean a fiat LCDM model 
with Quo — 0-3 and dark energy in the form of a cosmologi- 
cal constant, wxo = — 1. 
At the transition redshift 



~l.r 



-(l+3u*o) 



^MO J 



-l/(3tuxo) 



(18) 



the acceleration cancels out, R(z tT ) = 0. Then, we will fix 
the orbital parameters to their unperturbed values at z = z tr . 
In Fig. [T]we plot the time evolution of the orbital radius for 
an Earth-like planet (orbital period P ~ 1 year). Relative 
variations around our present time are really small. 

The behavior of an orbit depends on which energy-matter 
component of the universe stands out. Dark matter has an at- 
tractive effect and its density falls with time for an expanding 
universe. So, for a CDM model without dark energy, the or- 
bital radius expands as far as the scale factor expands. On the 
other hand, dark energy has a repulsive gravitational effect 
that can balance the dark matter action. Let us consider the 
condition for the expansion (contraction) of a gravitationally 
bound system, d(R/R)/dt > (<)0. The rate of fractional 
change of the orbital radius can be expressed as a function of 
redshift and cosmological parameters as 



x at 2 \ no I 



(19) 



x {(i + wxo)(i + 3w xo )n xl) R- 3w ™ + n M0 } . 

The fate of the bound system depends on the equation of state 
of the dark energy. In fact, if any dark energy is present, it 
will determine the expansion of the universe at future times, 
with no regard to the matter content. For either wxo > — 1/3, 
which means a decelerated expansion of the scale factor of the 
universe, or wxo < — 1, the orbital radius will expand with 
no regard to the total amount of dark matter and dark energy, 
see Fig. 12 For wxo > —1/3, the total source of gravitational 
attraction p+3p is positive and the dark energy is not repulsive 
anymore. The perturbation to the radius, which is oc R/R, 
stays negative and decreases in absolute value towards zero as 
the density decreases with time. The rate of change is as well 
decreasing, dx/dt — > for R — > oo. For — 1 < wxo < —1/3, 
dark energy is effectively repulsive (R/R > at late times). 
Its density and its repulsive effect drop with time. Then, the 
perturbation is positive and the orbital radius evolves towards 
its unperturbed value when the scale factor R diverges. In the 
case of dark energy in the form of a cosmological constant, 
wxo = — 1) the universe is destined to a de Sitter inflationary 
expansion during which both the orbital radius and its positive 
perturbation remain asymptotically constant. 

In phantom cosmologies, wxo < — 1, the dark energy den- 
sity grows with time and will eventually diverge in a finite 
time. The acceleration of the scale factor increases towards 
a singularity which is characterized by the divergences of the 
scale factor, the Hubble parameter H and its derivative H and 
the scalar curvature B17I1 . The repulsive gravitational action of 
the phantom energy will overcome the normal binding force 
of structures which will get dissociated lfl6l [l7ll . As for the 



case —1/3 < wxo < — 1, the perturbation is positive at late 
time, but both the perturbation, R/R, and the rate of change 
of the perturbation, d(R/R)/dt, grow with time. The orbital 
radius becomes larger with time with no regard to the matter 
content and it will get eventually unbound in a finite time. 

While the universe undergoes from a dark matter to a dark 
energy dominion, the orbit of a gravitationally bound system 
will move from a contracting to an expanding phase. The in- 
version redshift is 



-l + (l + u;xor 1/(3l0xo) (^r + l). 



(20) 



Such an inversion trails the transition from a dark matter to 
a dark energy dominated universe, z tr . As we have seen be- 
fore, the inversion would happen only for dark energy with 
equation of state —1 < wxo < —1/3. 

Let us now discuss the situation at the present time. The 
orbital radius is growing just now if 




<0; 



(21) 



the condition for expansion can be expressed in terms of cos- 
mological parameters as, 



^mo > -(1 + w X o)(l + 3w X o)^xo- 



(22) 



The equality between the sides of Eq. d22l guarantees for a 
constant orbital radius. We can easily see how, for our ref- 
erence LCDM model, a planetary system is expected to be 
expanding. This condition would change for JImo IS 0.2 and 



WXO 



> 



-1, see Fig. |2 



In the phantom scenario, the orbit radius will eventually di- 
verge. Then, the perturbation approach is not enough to track 
the complete evolution of the system, but alternative analyti- 
cal methods can be used ITol Flol . The dynamics of the radial 
motion of the bound system can be seen as determined by an 
effective potential. Assuming circular orbits, 



V cS 



-Xq 



(x/x ) 2(x/x ) 2 



(23) 

The equation for the minimum of the potential is the same of 
Eq. ( fLSb . The minimum with respect to x at time t gives the 
approximate orbital radius. The minimum exists only for lfl6ll 



J_i? 27 

nl~R [t > ~ 256' 



(24) 



with the equality picking out the time when the minimum dis- 
appears and the system becomes unbound. The dissociation 
time depends nearly exclusively on the late time behavior of 
the dark energy density. Then, we can find an approximate so- 
lution for the redshift of dissociation by considering £!mo = 0. 
We get 
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3wxo)^xo] 1/3 



(25) 
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FIG. 4: Rate of variation of the orbital radius in the Qmo ~ wxo 
plane in units of the rate occurring for a standard flat model with 
cosmological constant and Qmo = 0.3. We assume flat models. 



The span of time between today and when the system gets 
dissociated is 







dz 



(l + z)H(z) 



(26) 



As far as uq <C Hq, the dissociation redshift Zdis will be 
nearly indistinguishable from —1, with no regard to the prop- 
erties of the dark energy. Then, Afdis will be nearly indepen- 
dent from the orbital period of the system and the dissociation 
will take place nearly at the same time of the future singular- 
ity. In Fig. |3] we plot Atdis as a function of the equation of 
state. We can see that dissociation due to dark energy is not 
an impending danger. 



IV. ELEMENT PERTURBATIONS 

A perturbed planet orbit can be described through the or- 
bital elements of its osculating ellipse lfl8ll . The Lagrange's 
planetary equations account for time variation of these ele- 
ments. For a radial perturbation A r , we have II 1811 



da 

~dt 

de 
~dt 
duj 
~dt 



zA r sin if, 



e 2 A r sin if, 



i VT 



-Ar COS If, 



(27) 
(28) 
(29) 



where a is the semi-major axis, e the eccentricity, n = 
y/GM/a 3 the mean motion, lo the argument of periastron and 



if is the angle in the orbit plane from the periastron. Since the 
perturbation is radial, the inclination i and the longitude of the 
ascending node £1 do not change. 

In order to reveal the slow secular variation of the orbital el- 
ements due to the cosmological expansion, the above relations 
can be averaged over the unperturbed orbital motion, 



df 
~dt 



q(l-e 2 ) 
1 + e cos f ' 
Lo 



(30) 
(31) 



with Lq = GMa(\ — e 2 ). The average of a generic quan- 



tity / over an unperturbed period P 
through the integral 



27r jn is performed 



(1-e 2 ) 
2tt 



2^,3/2 rl-r; 



(1 + e cos if) 2 



(32) 



The cosmological term R/R can be Taylor expanded for an 
orbit as 



R 
R 



ip=0 




(33) 



Averaging over a period, we obtain the rate of variation of the 
semi-major axis 12311 . 



/da\ 
\dt/ t 




l + (l-e 2 ) 



( (1 



VvT 



(34) 

We turned to the time derivative of R/R thanks to Eq. OTb . 
The dependence on the eccentricity is pretty small, with a 
maximum decrement of ~ 25% for e ~ 0.3 and a null varia- 
tion from the circular case for very large eccentricity, e — ► 1. 
For a small eccentricity, 




|l-2e + 0(e 2 )} 



(35) 



For a circular orbit, e = 0, Eq. ( [34-b agrees with the result in 
Eq. ( fT3l . The zero order term in the expansion in Eq. (l33b does 
not affect a, which changes only for a scale factor with a non 
constant acceleration. In a model of universe dominated by a 
cosmological constant (de Sitter universe), the orbital radius 
does not change in time. This result can be also obtained in 
the framework of the Schwarzschild-de Sitter metric 1 1911 . 

Actually, variations in the orbital radius are pretty small. 
Over the life span of the solar system (~ 10 17 s) and assum- 
ing an Einstein-de Sitter universe with h = 0.5, the fractional 
change of an Earth-like orbit (P ~ 1 year) is of order of 10~ 23 
J5[]. Considering a span of time suitable for observations, 
~ 10 years, the orbital axis of a Pluto-like planet (M ~ M , 
a - 30 AU) changes by - 2 x 10" 14 m (~ 5 x 10" 15 m) for 
an Einstein-de Sitter model (reference ACDM model) with 
h = 0.7. The accuracy required to detect such small changes 
is well above the reach of next generation technological facili- 
ties. The today best accuracy on the orbital semi-major axis in 
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the Solar system, achieved for the inner planets, is of order of 
1CP 1 m 12011 . Variations get bigger for phantom cosmologies, 
but not enough. For a fiat models of universe with J7mo 5 0.1 
and wxo % —3, changes are 50 times bigger, see Fig. [4] 

Variations with time of the orbital radius are routinely ob- 
served in binary pulsars. As well known, orbital periods de- 
crease in time due to the loss of energy by gravitational ra- 
diation emission. Cosmology effects produce another mecha- 
nism that can affect the orbital period variation, 



P 
P 



3 a 
2a' 



(36) 



Let us consider the prototype binary pulsar, B 19 13+ 16 (P ~ 
7.8 hours). Its orbital period is observed to be decreasing with 
a rate of P = —2.4 x 10~ 12 with an uncertainty of < 10~ 15 
l2lll . In a flat model of universe with SImo ~ 0.3 and phan- 
tom energy with wxo = — 2, the orbital period variation is 
expected to be ~ 10 °/i 3 , well below the experimental accu- 
racy. The measurement of this effect in the near future seems 
very unlikely. 

The variation in the eccentricity can be expressed as 



1 VT~^ /da\ 

: \~dt/ 



2a 



(37) 



An expanding orbit turns more eccentric with time, as ex- 
pected from angular momentum conservation. As for the 
semi-major axis, changes are really small. 

A more feasible approach to detect the effect of cosmolog- 
ical acceleration on planetary motion should be to investigate 
departures from the third Kepler's law. As can be seen from 
Eq. ( TBI ), due to the perturbing acceleration term the observed 
angular frequency will differ form the mean motion n. It is 



5n 
n 



1 1 R 

~2n?R 
1 9 



(38) 
(39) 



where we have introduced the deceleration parameter, q, in 
the second line. At present time 



qo = i {0 M o + (1 + 3wxo)^Xo} • 



(40) 



The variation is now proportional to the acceleration of the 
scale factor. The statistical error on the mean motion for each 
major planet can be evaluated from the uncertainty on the 
semi-major axis, 5n = — (3/2)n5a/a, and can then be trans- 
lated into an uncertainty on the effective acceleration. When 
attempting to detect exotic physics in the Solar system with 
today data on changes in the mean motion, best bounds come 
from Earth and Mars Ull [l3ll . In our reference LCDM model, 
the today deceleration parameter is qo = —0.55. Consider- 
ing the Earth orbit (5a ~ 0.1 m, a ~ 1 AU), we get a limit 
\ a o\ ^ 0.8 x 10 10 /h 2 , ten orders of magnitude larger than 
estimates from observational cosmology. Constraints could 
greatly improve by considering a similar uncertainty on the 
measurement of the semi-major axis of outer planets, as could 



be obtained by dedicated radio ranging observations. For the 
Neptune orbit (a ~ 30 AU), an uncertainty 5a ~ 0.1m would 
imply |qo| IS 0.6 x 10 5 /h 2 , a steady improvement of nearly 
five orders of magnitude. 

The effect of cosmology on the periastron precession is of 
the same order of the change in the mean motion. All cosmo- 
logical models with a non null acceleration induce a variation 
in u>, 



(IlO 



R 



dt 2n R 



(41) 



In fact, differently from axis radius and eccentricity, the main 
contribution comes from the factor R/R and not from its vari- 
ation in time. At present time, in terms of cosmological pa- 
rameters, 



Tit 



2 H 2 q 



(42) 



In the case of a universe filled with only pressure-less 
dark matter with density pm and no dark energy, R/R = 
— 4ttGpm/3 and we retrieve the result in fl3ll . In a de Sit- 
ter model (Hmq = 0), then R/R = c 2 L/3, and we retrieve 
the well known result for extra-precession due to a cosmo- 
logical constant as usually obtained in the framework of the 
Schwarzschild-de Sitter metric Ull, [l9ll . 

The rate of precession predicted for the Earth is ~ —5 x 
10 _16 go^ 2 arcsec per year, ~ 10~ 16 arcsec per year in our 
reference ACDM model with h = 0.7. These estimates are 
ten orders of magnitude smaller than the state of art accuracy 
of ~ 10~ 6 arcsec per year obtained for the Earth and Mars 
orbits 1 22], so that the consequent upper bound on qo is of the 
same order of that derived from changes in the mean motion. 
We remark that dark matter and dark energy work in oppo- 
site directions. Constraints on the value of A from perihelion 
precession and mean motion in the Solar system have been ob- 
tained so far considering the contribution of the cosmological 
constant to the field equations but neglecting the contribution 
from diffuse dark matter 111 UN 2l and references therein]. The 
contribution of diffuse dark matter should be accounted for. 
Furthermore, the Solar system is supposed to be embedded in 
a Galactic dark halo with a density in excess of nearly five 
orders of magnitude with respect to the mean cosmological 
dark matter. This local enhancement makes the measurement 
of dark matter with accurate planetary astrometric data at the 
reach of future experiments lfl3ll . 



V. CONCLUSIONS 

Cosmological background can affect the orbit of a plan- 
etary system through the gravitational action of the smooth 
energy-matter components in which the system is embedded. 
Changes in the mean motion and anomalous periastron pre- 
cession are sensitive to the acceleration of the scale factor, 
R/R, whereas variations in the semi-major axis and eccen- 
tricity come from its time variation. In principle, a very de- 
tailed knowledge of the orbital parameters of a bound grav- 
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itational system would allow to put constraints on the cos- 
mological parameters through two diagnostics, i.e. qH 2 and 
d{qH 2 ) / dz, that are independent of the main ways in which 
observational cosmology usually investigate the energy con- 
tent of the universe, i.e. via either the distance-redshift re- 
lation or the growth of structure. This can be appealing on 
a theoretical side with its promise in breaking degeneracies 
which affect the other methods, but it seems an unlikely tool 
if we consider the observational prospects. Actual bounds on 
the deceleration parameter from accurate astrometric data of 
perihelion precession and changes in the third Kepler's law 
in the Solar system fall short of ten orders of magnitude with 
respect to estimates from observational cosmology. Even if 
we consider other tests of gravity physics, such as geode- 
tic precession and gravitational redshift, the situation should 
not improve substantially. In a previous paper 111 111 , we dis- 
cussed stellar system tests for the cosmological constant. Ei- 
ther future measurements of periastron advance in wide bi- 
nary pulsars or gravitational redshift of white dwarfs could 
provide bounds competitive with Earth and Mars data but this 
is not enough to constrain cosmological parameters. Better 
prospects come from future radio-ranging measurements of 
Solar system outer planets which could improve actual bounds 
by five orders of magnitude. Such an improvement can be sig- 
nificant but still not competitive with observational cosmol- 
ogy- 

Apart from the technological challenge required to reach 
the experimental accuracy needed to study cosmology in the 
Solar system, two other issues should be considered. In our 
analysis, we have considered a planetary system in an other- 
wise strictly homogeneous and isotropic universe. As a mat- 
ter of fact, planetary systems are embedded in galactic dark 
matter haloes whose density is in excess of several orders of 
magnitude with respect to the mean cosmological dark mat- 
ter. Furthermore, we have considered an homogeneous back- 



ground but tidal forces from the environment could easily play 
a role. These two effects on the galactic scale should easily 
overcome the effects of the acceleration of the cosmological 
scale factor. Any future observed deviation in the measured 
values of perihelion precession or in the third Kepler's law 
should then be seen as an indication for Galactic dark matter 
rather than an effect of an accelerated expanding universe. 

Till now we have considered effects due to the expansion 
only on planetary orbital motions. It might be interesting to 
consider also possible space-based interferometric type ex- 
periments similar to the one proposed to detect gravitational 
waves as is the case of the Laser Interferometer Space An- 
tenna (LISA) ll24ll . The radius of an Earth-like orbit would in- 
crease, assuming our reference ACDM universe with h = 0.7, 
by ~ 6 x 10~ 22 m per year. Accordingly, for a space- 
craft configuration like the one envisaged for LISA, which 
orbits around the Sun at the same distance as the Earth and 
with a separation of ~ 5 x 10 6 km between two modules, 
it would imply an increase of the separation of the order of 
~2x 10~ 23 m per year. In LISA a laser beam sent by a mod- 
ule bounces back from the other module's free falling mirror 
and then gets back to the original module where, after hav- 
ing bounced off the mirror, the light is mixed with a fraction 
of the outgoing light, and then interference is detected. In 
principle, an increase in the separation distance between two 
modules could be detected in this way but would require that 
the laser beam maintains a very high phase stability during 
years, which is well beyond today technological feasibility. 
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